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axis  of    no    feathering,    positive   aft.  Ibst 

I  i  M^ss  ii'  ment   of    inert  I  a  of   rotoi       I  i  I  out 

f  I  ,pp]  ng   hinge  .  lug   ft  ' 

Kj  Non-dJmen    lonal Izlng  term  Kj    =  pirR2(OR>2  I       . 

M  Longltudi     1 1    p  i  tch  i  ng  m<  ■  -cut   the   c .    . 

of    the   aircraft,    positive    nose    jp.  ft    I 

Mf  L  ngitudinal    pitch]     ;  morm  n1  to  fus    I 

^  and   powt  r  effects  ft    I 

My ,   Ma  1 1  I  zed   pit    h  1  ng  mom  >nt    der  I  v   1  i  . 

Np  Normal    force   due  to   thrustli    |   propeller, 

positive  up  lbs. 

0  haft    torque 

Main    rotor   radii  ft. 
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ti     _■   j-i      Vcosa 
u,  Tip  speed  ratio,  u-  = — ^-  n.d. 

p  Mass  density  of  the  air  slugs/ft 

a  Rotor  sol Idlty  a  »  ^£  n.d. 

ttR 

fi  Rotor  angular  velocity  rid/sec 

To  indicate  that  rotor  speed  is  a  varjal  |e  In  accordance  with  torque 
equilibrium  derivatives  so  taken  are  designated  with  I  _  ~o  Derivatives 

ly^U 

taken  with  rotor  speed  specifically  held  constant  will  be  de  I  ,n  ted 
with  |  . 
IP 
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SUMM  RY 

In  an   I  mpt  to  quantitatively  exami  ie  1   effect  of  variable 
rotor  speed  on  hhe  loi  jjtudlnal  static  ;1  i  Ility  of  in    horot  tl 
rotor,  express  ions  have  von    developed  for  the  changes  h  rotor  forces 
under  a  condition  of  mafntei     of  const  nt  zero  torque.  R 
moment  derivatives  were  then  calcul  ited  for  i  small  >utogyro  rotor  u 
two  limiting  conditions,  l«e.  for  rotor  speed  assumed  c  i  I  n1  I     ih- 
out  velocity  and  attitude  pertui  >ati.ons  and  for  roior  speed  v  ;  ry  i  g  In 
such  a  way  that  equilibrium  torque  is  maintained  throughout  the 
perturbat ions,  Comparison  of  the  two  conditions  showed  thai  illowi  ig 
rotor  r.p.iru  to  vary  caused  the  velocity  stability,  My,  to  go  from  a 
strong  stable  condition  at  constant  r.p«m.  to  a  neutrally  stable 
condition.   The  effect  on  .ingle  of  attack  stability  was  to  cause  an 
unstable  condition  at  constant  r„p.rr, .  to  become  stable* 

.■tic  stability  derivatives  were  calculated  for  a  small  iuto- 
gyro  under  the  variable  rotor  speed  conditions  and  steady  state  flight 
tests  were  conducted  to  attempt  to  verify  these  derivatives  from  trim 
curves.   Although  the  flight  tests  tended  to  confirm  the  results,  they 
were  felt  to  he  inconclusive  ]n  that  they  could  only  verify  ratios  of 
the  derivatives*   Using  the  const  >n1  coefficient  of  the  char  cteri  Mc 
on  "frequency"  equation  as  a  measure  of  static  stability  the  ilrcraft 
was  determined  to  be  statically  s1  ible.   Further,  the  v  ;r  i  itjon  of 
rotor  speed  improved  the  static  stability  although  It  did  reduce  the 
ve loc  i  ty  stabi I i  ty. 
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It  was  further  concluded  that  the  static  stability  of 
the  aircraft  was  independent  of  horizontal  e.g.  position  j,  long 
as  the  range  of  e.g.  movement  was  confined  to  control  ||mI1 


AN    INVESTIGATION  OF    THE    LONGITUDINAL     TTIC    STABILITY 

OF   AN  AUTOGi 


I.   INTRODUCTION 

It  has  been  stated  and  is  generally  known  that  the  artli  u  I  ited, 
lifting  helicopter  rotor  in  forward  flight  is  statically  stable  with 
respect  to  velocity  changes  and  is  statically  unstable  with  re  pecf  to 
change  of  angle  of  attacks  Th ]        : 1  implies  the  cl  issical  rigid 
rotor  witl  no  artificial  type  of  stabilization  devjo      fui  rh  i  is  wmes 
that  the  roi  itjon  il  spec  J  of  the  rotor  is  constant  <, 

*   ssow  and  Myers  (Ref.  I)  state  hh  I  the  lutorotatlng  rotor  has 
different  stability  characteristics  than  the  powered  rotoro   The  primary 
reason  for  this  difference  Is  the  fact  that  the  rotor  speed  of  the 
autorotating  rotor  is  not  controlled  by  the  engine,   ut  Is  free  to  vary 
with  *  h  r,  ]c  In  forwai  d  speed  or  angle  cf  ittack.  Hohsnemser  (Refo  2) 
states  1  hat  the  effect  of  these  vari  tions  In  rotor  speed  is  to  make  the 
autorotating  rotor  neutrally  stable  with  ch     i  in  speed  at  constant 
angle  of  attack  and  positively  i1  il  le  with  change  In  ingle  of  attack 
at  co.'i '.tart  speed. 

It  is,  therefore,  the  objective  of  this  paper  to  derive 
expi      s  for  the  longitudinal  static  stability  of  the  autorotating 
rotor  in  forward  flight,  to  examine  the  effect  of  allowing  the  rotor 


sp  ed  to  v  iry,  md  to  ittempl  to  v  rify  \   e  r<  ults  111  Hvely, 

If  nci  quant  11  't  i  ^c  ly,   y  s  imp  |e  equl  H  trium  f  1 1  ;  I"  Ing  of  i 
small,  one-     ut  ^yro.   See  Flgun  l<> 

Throughout  the  analysis  stability  derivatives  will  be  kept 

in  ci  jlmenslonal  fen:.'  to  retain  the  physical  conce'pl  ;  md  the  phy  i 

quantities  of  the  test  aircraft,  as  li  1  I  i  i  Appendix  \,    will 
he  used. 


II.   MALY3  I  S 

In  the  stability  analysis  of  ar y  vehicle,  it  Is  first  necessary 
to  establish  a  system  of  coordinates  to  describe  the  notions  of  the 
vehicle.  Since  actual  flight  testing  Is  involved  In  this  analysl  ;, 
it  is  felt  that  the  "stabl  I  Ity-axl  i"  y  tem  where  the  origin  vf  coor- 
dinates is  placed  at  the  e.g.  of  the  aircraft  and  the  po  Itjve  X  axis 
fixed  in  the  aircrafi'  and  initially  oriented  into  the  re  I  at  i  ve  wind  Is 
best  suited.   Flgure'2  shows  the  axis  system  and  various  quantities 
relative  to  it. 

How  to  accurately  describe  the  motion  of  the  vehicle  depe 
upon  the  number  of  degrees  of  freedom  it  has  ard  how  strongly  each 
degree  affects  the  motions.  As  it  Is  our  purpose  to  examine  the  longi- 
tudinal static  stability  and  to  determine  the  effect  of  the  rotor's 
rotational  degree  of  freedom  upon  It,  certain  simpllfyi      jnptions 
re  necessary  lo  avoid  making  the  calculations  so  complex  that  an  engineer- 
ing analysis  becomes  Impractical.   It  Is  first  assumed  that  the  longjtudlna 
and  lateral  motions  of  the  aircraft  are  loosely  coupled,  so  that  the 
longitudinal  motions  may  be  effectively  described  by  the  classic  three 
degree  of  freedom  equations,  I.e.  X  Force,  Z  Force,  and  Pitrhing  Moment. 
The  rotor  itself  has  six  degrees  of  freedom,  I.e.   coning,  Literal  and 
longitudinal  flapping  and  steady,  longitudinal,  and  lateral  twlstln  , 
plus  the  freedom  of  rotational  speed.   Under  the  assumption  of  3  rigid 
rotor,  which  is  usually  valid  for  stability  analysis,  and  writing  the 
rotor  statics  In  terms  of  fuse  I  a je  variables,  the  rotor  degrees  of 
freedom  may  be  reduced  to  the  variation  of  rotor  speed. 
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Further,    sinc<       t<    \d)       I     I       I  I  1  ih !"    I  hn  Jques      r 

comp  I    I     ly      ultjble    for  tl  I       '      !       Hve    ,    rate     ind  ao     I    i    itloi 

1  erm     will      e  <  m  Itted   from   1  |ujt  Ion  Hon. 

Z   Foi  uatloi 

lummjng    forces   along  the  Z   jxIs    (see    Figure      )  glectji 

any    fu  ;e!age    lifi    as   small    compared    ho    I  he   roti  r    for   es    gl 

ZF     =  .Vcos-y  -  Tcosa  +  Hsjna  (I) 

F        sma  I  i      n    i     \}    the   cosine  may   be   t  ikon   as  equal    to   u  iltyo    II  i 

!!    I    i    i!     I"  he  |  r  of    small    quant  I  tie    ,    the  Z    force   equ   tlon   reduces 

to 

EF7   =  W-T  (2) 

Tak  ing    i  i    r   m      ra Is 

ZF.   +  1F2    =  W-(T+  T)  (    ) 

5ul  tlon    (2)    f n  'i    (    ) 

SAF2    -  -AT  (4) 

In    i erms  of      he    i otor 


T               T  T                   T                ?>T 

T     =  —  v   +— (Acli    +  A6)    +t a,    +  —  &-. .    + :    , 

V  oa         |              ,       I          I    .       1 

T      •  5T              T 

,     '  da                                           (5) 


he  approximation  i   mad'  th  I  the  rotor  Is  restrained  to  neve  in 
he  pi  nie  of  symmetry  In  determination  of  rotor  thi   I     that  the 


. 


I  vat  1  ve     . ,      i  ■:.   i  I     '        ma  I  I     is    is       iggi  I    ii  I  .      }    1 


i     I     I  .  !■"   p  ' 


^T      ?J  r T 

AT  =  — AV  +  — (Aa.  +  A6)  +  — ..0 

•V  I 


(6) 


hi 


|n<  I  led  In  equation  (6)  Is  the  term  ^qAP   Implying  the         f  n  hor 

speed  as  an  independent  variable.   Examination  of  ;he  torque  >    Non  for 
an  au k rot  t In  ;  rotor 


blrotor  -  +  Qaerodynamlc  =  0 


(7) 


shows  that  for  steady  state  conditions  the  aerodynamic  torque  mu  \   equal 
zero.   If  the  assumption  is  made  that  the  Inertia  of  the  rotor  is  negligible 
then  the  speed  variation  of  the  roior  is  dependent  on  torq     illlbrlum» 
This  assumption  allows  us  to  make  rotor  speed  a  depei,de:  i  variable  and 


eliminate  the  torque  equation  from  the  equations  of  mctlc;u   Fur1   i 


defining  the  rotor  force  derivalives  in  a  special  way  such  fhat 


dT 
3a 


Q=0 


means  the  rate  of  change  of  rotor  thrust  with  angle  of  „ittjck  while 
forward  velocMy  is  held  oconstunt  but  roior  speed  Is  allowed  to  v_iry 
In  accordance  with  torque  equilibrium,  then  equation  (6)  may  be  wrI1 


.  ii 

T 

T 

AV   +    — 

=   AV 

9=0            S<X 

Q=0 


(Aa,  +  A6) 


(3) 


and  the  summation  of  force  changes  In  the  Z  direction  for  equilibrium 
f I i  ght  i  s 


. 


SAF- 




iV 


0=0 


AV  + 


dZ 


(   ,  +  A6) 
Q=0   ' 


(9) 


V 


?-0 


«v 


g=o 


(9  ) 


_ 
da 


g=o 


dT 


0=0 


Pitching  Momeni  Equation 

iln  referring  to  Fig.  2,  the  moments  about  I     ircraft  c«g« 
may  be  written  as  follows 

EMcg  =  x(Tcos6  -  Hsln6)  -z(Tsln£  +  Hcos8)  +  Mc;  +  Mf     (10) 

Since  The  rolor  on  l"he  alrcrafi  being  |i  .   ,I  jated  is  a  s  -saw  typ 

the  momeni  due  to  flapping  hinge  eccentricity,  Mc,  Is  zero.   This  moment 

fherwise  often  be  neglected  if  the  f  lapping  hinge  offset  Is  r.i  II. 
Further  assuming  small  angles  and  neglecting  H6    1  e  product  of 
titles 


EMcg  =  xT  -  z  (T5  +  H)  +  Mf 


(II) 


Tak  in  i  j  nc rcne'   I  )  I 


EMcg  +  M, 


g  -  x(T  +  AT)  -  7  [<T+.T)(6+A6)+<H+AH)]  +  Mf  +  Mf 

P     'p 


(12) 


lubtractlng  equation  (||)  from  (12) 


SAMcg  =  xlT  -  2  (ATS+A6T  +  J\)    +  JAf 


(I  ) 


i      (t  Itutlng    i  i  ro  equ  i1  ion    (  I   0 


•T   =SLL 


•I 
V 


V    +  SI I        (M,    *  A6) 
x,=o  5a  Q=0 


(  ) 


!! 


11 
Mf 


Q-0 


.V    +  _ti  <£K.    +  Li) 


5Mf 


_M. 


v  +  — _Aa, 
1 


(i     ) 


(I    .  ) 


yle  I  ds   the    f o  I  low  I  ng 


5T 
-"Meg   =   (x  - 


5T 

-z5  — 

Q=0  ^V 


5H 


-z 


5T;  ST 

+  (x  —         -z6  — 

Q-0  da 


Q=0        jV:Q=0        oV 

5H  I  5Mf 

-2   — !  + 

Q=0 


B.  )  v 


-z  — I  +  

da rwi       Sa 


P     1A 


) 


jQ=0 


(14) 


da 


-zS 


Q=0 


5  a. 


-Z 


H 


Q*0 


a 


00. 


-zT) 


Q=0 


Considering   equation    (14) 


5M 


■eg 


njV 


AV   + 


Q=0 


5M 
5a 


Q=0 


Aa,   +  — 
56 


A6 


Q=0 


(15) 


where 


5M 

5V 


Q«0 


5M 

5a! 


Q*=0 


5M 
5V 

5M 

oV 

5M 
da 

5M 
5a 


+  ^k 


Q=0  |  rotor 


5  V 


.5T;  T 

Q*=0      rotor  v    Q=0  v 


-zill 
Q=0     ^V 


Q=0 


5M 


Q=0      rotor 


3a, 


Q=0;    rotor         JU| 


-  2l 
5a 


5T 


-z 

9-0 ' 


-z- 


Q=0 


5H 

"5a 


(  !      ) 


( ;     ) 


(  ! 


(I     :) 


s 

'  0=0 
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-  .1 


(I  ) 


0=0  I  ( ( 
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I  i    of  Mot  Ion 

Tii1  previously  developed  equations  of  motl<  n  I  r   I      hate 
fli  :,i   onditions  can  row  be  written 


V 

v+_ 

.v;=0 

_11 

M 

V    +  — 

V  ;  ^ 


r'    i; 


0=0  |Q«0 


+  ^** 


=  0 


q=o 


Ac  =  0 


( 1 6a ) 


(I   ) 


0=0 


And  H  Is  possible  to  make  a  qi  illtatlve  check  on  fhe  derivatives  by 
comparing  the  slope  of  the  trim  curves  taken  from  steady  state  flight 
test  to  the  predicted  slopes^  si 'icy  they  may  now  be  predicted  from 
equal  ions  16. 
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A  check  on  the  rotor  thrust  derivatives  may  i|so  be  made  l>y  noting 
from  equation  (I6a)  (since  t£L  -   Adj  +  AS) 
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Further,  the  effect  of  variation  cf  rotor  speed  on  each  partlcul  -r 
der  ivative  is  available  by  computing  each  one  separately  under  the 


assumption  of  constant  rotor  speed,  letting  th«  rotoi 


derivatives  — 
iV 


a 


■a 


=  0,    and   drop    those     I 


Evaluation   of      he  Derivatives 

From  the   definition  of    the   rotor   thrust   coefficlen 


T  «  PTTR2(nR)2  Cj 


(I    ) 


where  C  is  the  trim  value.   Letting  Kj 
1  V 


=  pTTR2(GR)2,  then 
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Including  Cl   as  an  independent  variable  reqjjres  the  addition  of  the  torque 
equation  to  the  equations  of  met  ion.  As  we  have  already  seen,  however, 
neglecting  the  inertia  of  Jhe  rotor  allows  us  to  make  0  a  dependent 


variable,  and  the  rotor  thrust  derivatives  could  be  expressed  as 
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(201  ) 


Ev  lu.tlon  of  the  thrust  coefficient  derlv.tlves  may  be  made  fro;   i 

element  theory,  where  the  thrust  coefficient  is  giver;  by 

r.2   .,2 


30  13     2  J       2    4 


(2  1) 


10. 


For  a  constant  blade  pitch  angle  the  thrust  coefflcienl  Is  i  function 

of  the  variables  p,  and  A 

Cj  -  Cy(u.,X) 

f 

If  the  usual  assumption  for  the  average  Induced  velocity  cf  a  lightly 
loaded  rotor  In  forward  flight  is  taken  from  momentum  theory  as 


v  =  ! (22) 

2prrR2V' 


where  V'  Is  the  resultant  velocity  at  the  rotor  disk,  then  the  inflow 
ratio,  A,  m>jy  be  expressed  as 

CT 
2  7^*2 

Since  by  definition  the  tip  speed  ratio,  u,,  Is  given 

Vcosa  .->a\ 

p.  ■  — — —  (24) 

i  R 

It  then  follows  that  the  thrust  coefficient  could  be  expressed  as  a 
function  of  the  variables  V,  a,  and  ft. 

CT  =  CT(V,a,n)  (25) 

The  thrust  coefficient-  derivatives  may  !hen  be  evaluated  following  the 
chain  rule  for  partial  differentiation 

5£i  =  ^£T^  +  2El!i  (26  > 

>V     u.  v   Sa  3V 
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Note   th  to   the   rol  ition    between    inflow   ra1  io     n  d   thru    t 

coeffi    i       I    fhrougl    the  momentum  theory,    l"he    Inflow   ratio    jerlv    M\ 
are   a    little   compile  d   have   a   ^h   r.     er   somewhat   more    like  icti 
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12. 


L>ome  simplification  i:  pc  ■  i  l<  i  iii  1 :  f  low  i  itlc  lerlv  itlv 
:  y  ne  )l ec1  Ing  \z  comp  ired  to  pr  at  1  ip  p  I  r  [  3V  0»  I  i 
and  approximating  the  inflow  ratio  as 

CT 
•  =  M.0  -  — 

Sir, ce  much  of  the  range  of  flight  being  investigated  vas  below  p  =  0«  I  , 
thi>  simplification  was  not  used  in  the  included  calculations. 

Evaluation  of  the  rotor  rpeed  derlv  fives  musl  evolve  from  tl 
torque  relation,  since  it  Is  the  change  of  rolor  speed  in  accoi  I     .vith 
torque  equilibrium  1h=t  is  desired.   From  B  I  ley's  aerodynamic  torque  re- 
lations (Ref.  4)  the  aerodynamic  torque  m  .y  e  expressed  as  a  fu-ctlon  of 
the  tip  speed  r^iio  jnd  inflow  r^tioo   For  a  fix!  blade  pitch  ingle 
then : 

CQ  =  C>/V>'-) 

See  Appendix  II  for  complete  torque  relations.  Since  the  only  ton 
acting  on  the  autorotatjng  rotor  ;t  equili.rium  Is  aerodynamic  torq 
it  f o I  lows  that 

Cq  =  0 

and  the  inflow  ratio  may  ne  determined  In  rerms  of  the  tip  speed  r^tio. 
Tills  ha;    i  n  comp  ited  for  trie  rotor  being  Invi   i  |  ted  in  \    Is  plotted 
on  Fi  gure  3. 

From  the  definitions  of  the  tip  speed  ratio  ii  I  h  flow  ritio  the 
aerodynamic  torque  coefficient  may  :e  expressed  as  a  function  of  the 
v.ir  lab  |es  V,  a,  and  0. 


I  . 


I i  fol lows  that 
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v    iV     da     on 


'An 


H  v  i  j  neglected  the  Inertia  of  the  rotor  is  h:vi:  j  •  negligible  effed 
on  rotor  speed  perturbations,  The  equlli  t-rlum  orque  mus1    maintained. 

Tli  on 

ACq  *   0 


so  ; hat 
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oCq 
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5n" 


These  particular  derivatives  may  be  evaluated  by  direct  differentiation; 
they  have  been  so  performed  for  the  rotor  being  Investigated,  and  aru 
plotted  on  Figures  4  and  5  as  functions  of  the  tip  speed  r  i1  lo,  p,. 

Substituting  equations  (27)  ant  (2j)  back  into  (20)  determines  the 
rotor  thrust  derivatives  under  the  varying  rotor  speed  conditions*  These 
derivatives,  alony  with  thrust  derivatives  computed  on  the  assumption  of 
constant  rotor  speed,  have  been  computed  for  i he  rotor  being  investigated 
and  are  plotted  on  Figures  6  and  7. 

Since  the  trim  values  of  the  inflow  rcitlo  X   can  be  calculated 
from  the  aerodynamic  torque  relation  In  terms  of  the  tip  speed  r<j"rio  u,, 
as  Is  shown  in  Appendix  II  and  plotted  on  Figure  3,  the  trim  thrust 
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coe-f  M<  i  •■  nt,   C-p    i  rm]i      I    from     i  I     I        lei        I     fheory,       |uat  loi     (2I)« 

Th  i  n     loi  e  and  the  'trim     hrust     oef  f  i  ciei       !  .    p  Ic 

ioi     >f   the  tip      peed   i     M  Figur     8.      For         liven  gross  we  I  , 

the   trim   re tor    sp       I   may     e   ob1     i ned    fron    >  q  i   t I o      (I    ),      The    trim 
v  i  I  ue  of   1  he  angl      c  f  i    fh<     p  Jane      f   no   feather  I         i       I        r- 

mlned    from  equatloi     (2>),     md      rim    f o r vv  \\  d    speed    is  '  r  he 

deflriti  in  of    l"hi     hip      peed    ratio,    equ      loi     (24). 

Evaluation   of    th-   momenl     lerlv    five    ,    equation      (I    ) ,    requli 
a   knowledgs    of    fhe   II    force   d<  rlv    five    ,    trim   control     in  \\      ' ,      nd 

lage   and   powf-r   contributions  I  I  I  ^  V  ^    AMf    >     ■■   well     is   the 

rotor     fhru    r     i    <r\  i   tiv  Th''   rotor  H    foi  oeffl    lent  jivei     by 

Gessaw  and  Myers    (      f.    I)    for    hhe  H    force    ii       he   pi  me  of    no    feath   rl 
Is 
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a  a 


M   =  _       ^\ 

'    la  5 


,  e 
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■     I  M-  '\  -o     |  M.--Q 

2  4  4  4 


(       ) 


Unce  the  fl  >ppi"j  angles  relatjvi  to  the  plane  of  no  featherli     n 
expressed  as  fui  Mons  of  the  tip  speed  ratio  and  inflow  ratio  (the  coning 

lie  for  the    -  aw  rotcr,  a  ,  Is  built  i        nstant)  they  may 
substituted  Into  the  H  force  coeffji  i     quatlon  and  it  may 

completely  expressed  as  a  function  of  p,  md  -. 

o         a 
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CH  ■  CH  (p.,  '-) 


I ;  v  i     I" a  1  ned  >      is  a  f  u net  J 01   f  p,  f  r  om  toi  |  i  equ  I  1 1  br  I um,  the 

trim  H  force  coefficien   Is  tl          i    i  ,            plotte 

in  Fi  ure   ..   Furthi  r,   knowji  ]  the  trim  r  tei  p   I    erml        I! 

force  Itself  from  fhe  definition  of  the  H  for  ffl  lent 
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V  pnR2(i  R)2 


(  0) 


The  H  force  his  been  calculated  for  the  rotor  under  |i     '  |atIon 
is  plotted  Dn  FI  jure  9  as  a  f unci  Ion  of  tip  speed  ratio  u* 

From  the  definition  of  the  M  force  coeffl  lent,  equation  (30), 
the  li  force  derivatives  may  be  obtained  in  a  slmll  jP  m  mi  er  to  thit 
used  for  the  thrust  force  coefficients. 
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Evaluation  of  the  H  force  derivatives  under  conditions  of  varl  il  le  rotor 
speed  also  follows  simil  rly: 
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10. 


Since  the  H  force  coefficient  i  n  <pi     i   completely  [n  ferms 

-  f  the  tip  speed  ratio  u,  and  inflow  ratio  \    by  substitution  of  the 
flapping  angle  relations,  equations  (29)  ^  into  equation  (21) }    it 
fo  I  lows  tli   : 

CH  =  CH(V,~,'  )  (33) 

The  H  force  coefficient  derivatives  m  y     'valuated  by  differ  ntl  tlon 

i         r  siml  I  jr  to  th  :t  used  to  obtain  the  thru  ,t  force-  coeffl  i-  t 
der  i  vat  i  ves, 

The  rotor  H  force  derivatives  have  been  calculated  for  both 
constant   r.p0m.  and  variable  r.p.^.  as  indicated  above^  and  re 
presented  on  Figures  10  and  IK 

To  find  the  trim  contra  I  angle  5  it  is  necessary  to  solve 
the  moment  equation,  equation   (II),  for  trim  conditions 

EMrc.  =  xT-z  (T6-+H)  +  Mf   -  0  (34) 

cy  p 

which  necessitates  an  estimation  of  the  moment  about  the  c«je  due  to 
fuselage  and  power  effects,  Mf  «   To  a  first  approximation  it  Is  f<  It 
that  the  fuselage  itself  will  have  little  aerodynamic  moment,  If  >ny« 
It  does,  however,  have  a  small,  flat  p  |.ite  horlzor.t-il  tall  surface 
which  may  offer  an  appreciable  moment  at  higher  angles  of  attacko   The 
tail  contribution  to  trim  moments  would  be 

M-f-  =  Xj-S-j-a^-Oj-q  (35) 

Power  effects  are  difficult  to  estimate  accurately.  A  small  amounl  Is 
contributed  due  to  direct  propeller  thrust  as  the  thrust  axis  is  offset 
a  small  distance  from  the  c.g0 
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MTp  -  ZpTp  (  '  ) 


A  propel  lei    rma  I  foi  e  will  exi  I"  wh  n  I"  he  thru  t  axis  (parallel  to 
3  fuselage  n  feren  e  line)  op<  rates  at  an  mgl  of  att  icko   Followli 

j  of  Rlbn  r  (Ref.  5)  i  propeller  normal  force   ised  on  prop  I  ler 
thrust  coefficU      1  pn  p>  I  ler  geometry  has  beer  i  -tlmated.  The  nor 
duo  to  propeller  normal  force  can  th<  n  be  expressed  is 

M'JP  =  XP   — "  >  (  7) 

Since  the  total  pitching  moment  due-  to  fusel  ige    :  power  effects  is 
now  ava  i  |a  k  _, 


M   =  M  |  +  Mr   +  M, .  (38) 

rP         P     P 


not  I ng  th a  r 

I  =  ~f  =  ~p  "  CC-6-lR  (39) 

the   trim   control     angle    6   m  ,y   now    be  obtained    .  y   su  tstltutjng   equation    (38) 
i  i  t      (  34  )    and      rear  ran  3I1  <ju 


xT-zH-f*+  ;t  i+q(a-JR)+zpTp+xp—  (cx-Ir) 

-  ^r 

Zi+X-j-  j-|.  j-j.q+xD Si 


(40) 


Note  that  the  power  moments  require  an  estimation  of  the  hhrustli 

propeller  force,  which  may  be  obtained  by  solving  the  horizontal  trim 
force  equat]   : 


vp     _  TpCosap-Tilna-+lcosa-Df-Ws!nY-0  (41) 


imi         si  i  1 1  !     •     lies 


Tp   -  Ta  +  H  +  Df   +  (42) 


f-jroatlon   of    t  i   trlbutlon   of    fuselagi     and   power   effects   to   the 

M      stability  of   the   aircraft,     Mf    ,    will    probably    Introduc       i      great 

P 

an  error  as  an\  othei  approximated  .  As  was  Indicated  In  the  analy  1  -  for 
the  trim  c   rol  angji  ,    ; ,  the  fuselage  Is  assumed  to  have  a  negligible 
momenl  of  Itself.   No  known  theory  is  avail  ible  to  predict  a  moment  for 

i   fr  me-lll   structure  and  without  wind  tunnel  data,  It  can  only  be 
assumed  small.  Tall  contri  utlon  to  stai  llity  can  be  estimated  from 
equat  loi  ( 3  )  : 


=  x-|-3j-  i^-a-f  V  (43a) 

— -  =  x+Sta+q  (43b) 

be. 


.  Ince  the  momenl  arm  of  "the  thrust  axis  Is  quite  small  and  no 
simple  theory  exists  for  estimation  of  f"he  engine  thrust  changes  with 
velocity  f  i  constant  pitch  propellers,  tlu  direct  propeller  thrust 
contribution  Is  also  neglected.   The  propeller  normal  force  contribution 
to  stability  is  "taken  from  equation  (37): 

cMNp       BNp 
=  x   


The  hall  ind  propeller  normal  force  contri  ut I  -   to   I   lllty  h 

Li  I  ited  as  Indli  ted  ibove  and  ore  presented  on  Figur   I  . 
,','ii!i  i      hes  of  The  rotor  force  dei     Ives,  trim  control 
le,  ii  d  fuselage  and  power  coi  tri  utlons  a.  i!  i  |e,  equatloi   (I  ) 
may  now  be  solved  for  the  moment  derlvatlveSo  The  J-  rlvatfves  hcive  Seen 
obtained  for  ;  he  rotor  alone  for  both  constant  rotor  speed  and  v  ri  i  I 
rotor  speed  conditions  In  accordance  with  torque  equilibrium      re 
presented  or    Fi  jures  13  and  14.   Further,  the  aircraft  derivatives  have 
been  calculated  I   I  <  ludlng  the  fuselage  and  power  effects  and  )re 
preserted  or;  Figure  i  15,  \y>,    and  17. 


III.      EQUIPMENT   AMD   PROCEDURE 

I 

Ai  rcraft 

The  aircraft  used  to  obtain  trim  •    I   for  thi  .  anaiysl  .   was  i 
sma I  1 ,  one-man  autogyro  known  as  a  Bensen  B-8M  "Gyrocoptero"  See 
Figure  I.   The  aircraft,  built  and  flown  by  the  author,  Is  powered  by 
a  72  hP  McCulloch  air  cooled,  four  cylinder,  two  cycle  en      riving 
a  fixed  pitch  Banks  and  Maxwell  pusher  propeller..   The  rotor  |s  a  see- 
or  "teetering"  rotor  with  a  diameter  of  20  feet  and  a  constant  i  lade 
chord  of  6-3/4  inches.  The  blades  are  untwisted,  have  a  fixed  blade 
pitch  angle  of  1 2  degrees,  and  a  built-in  coning  angle  of  4  degrees. 
Control  is  affected  throujh  an  overhead  control  slick  connected  directly 
to  the  rotor  shaft  by  two  bearings.   The  rotor  shift  I  ,  free  to  rotate 
10  degrees  laterally  and  I  or  1  Itudli  My  about  the  main  thru  t  bearing  I 
the  control  head  (rotor  reference  axis).   A  conventional  vertical  fin 

rudder  Is  Immersed  in  .he   propeller  slipstream  at  the  tall  to  con- 
trol sideslip.  A  sma i I ,  fixed,  flat  plate  horizontal  tall  surface  Is 
mounied  directly  under  rhe  engine  aid  functions  both  as  stabilizing 
surface  and  propeller  shield  when  raxlin  ,.   The  aircraft  w>  i 
225  pounds  empty  and  4C0  pounds  at  normal  gross  weight, 

I  nstrumentat Ion 

Ai  rspeed  and  A  1 1  i  tude 

A  standard  helicopter  airspeed  Indicator  and  sensitive  altimeter 
were  used  to  determine  airspeed  and  altitude.  The  I   tr    ■  Here 
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drive     )  swivel  I  In  |   I  lie  pressure  probe  am     hielded  toi  )|  bead 
pick-up,  both  mounted  on  a  small  boom  on  the  forward  portion  of  the 
aircraft.   The  airspeed  system  was  calibrated  by  the      lard  sp< 
course  method  and  the  calibration  curve  i  :.  pre    I  j  on  Figure  I  ■> 

Rotoi  '..peed 

tor  speed  was  obtained  by  driving  a  Servo-Tek  Pec 
linear  tachometei — generator  with  a  short,  flexible  ihaft  directly  from 
the  rotor  shaft  and  was  displayed  on  a  shunted  micro-ammeter  mou  I   on 
the  Instrument  panel.  The  system  wos  calibrate     driving  the  tacho- 
meter-generator with  a  variable  speed  drill  motor  and  measuring  r.p.  1. 
with  a  stroioscope.   The  micro-ammeter  was  then  biased  so  that  ir  ro^'i 
rotor  r.p.m.  directly.  The  accuracy  of  the  system  was  within  the  i  Ility 
of  the  pilot  to  read  the  ammeter,  which  is  felt  to  :e  +  2  rop.m.   Rotor 
r.p.m.  was  used  only  as  a  check  on  the  predicted  thrust  coefficient  . 

Control  Angle 

Control  angle,  corresponding  to  stick  position,  was  measured  by 
mounting  a  7 50  ohm  potent iomctcr on  the  mast  just  below  the  rotor  control 
head.  A  waxed  nylon  line  wrapped  several  turns  around  a  pulley  or  the 
potentiometer  shaft  had  one  end  connected  to  the  aircraft  frame  throujh 
a  tension  spring  and  the  other  end  connected  to  a  short,  rigid  arm  mounted 
on  the  movable  rotor  head.   The  tension  provided  by  the  spring  :||owed 
the  wiper  of  rhe  potentiometer  to  follow  the  longitudinal  movement  of  the 
rotor  shuf 1  „   The  potentiometer  was  driven  ly  two  1.34  vol  I  mercury 
.arteries  and  the  output  was  displayed  on  a  micro-ammeter  mounted  c 
the  instrument  panel.   The  system  was  calibrated  by  measuring  control 
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!  •   .  |  propcllerprotr  i  i       i         I  me 1 

i0,000  t  hm   pot      !  iomi       r    i  i   'S   wl1       It.      The 

ca  I  i  to  n  I  I  y    i       legrees, 

of    t"h<  i    m   ttas   within    the   pilot's   ability  to   read    N  r, 

is    i'   l  .     I        e  approximately   +  0. I     legre    < 

Fu  ,ngl  ?  of 

Angle  ol    I  i  k  •   1  l"he  fusi  I  i  ;■•  reference  line  (  ee  Figure  2) 
isured   y  a  small;  balanced  balsa  wood  v  in         on  the  for* 
extremity  of  the  ah  :raft.   The  vai  e  w  is  mou  h  I  on  i  horl  ontal  shaft 
rotatli  |  i   all  bearing  and  a  polntei   i   e  ted  to  the  sh.ift  r 
fuselage  angle  of   ft  »ck  or  a  scale  oi  ("he  vane  frame.  The  system  was 
call  rated  with  a  propeller  protract  r,     md  Its  iccuracy  Is  felt  to  be 
within  +  0.5    i  es. 

is  .■   i  ih1   .I  Cen  er  of  Gravity  Determination 
The  ro  s  w<  Ight  ind      r  of  n  ivjty  (e.g.)  of  the  ijrcraft 
was  i  :  ree  different  methods  as  a  cross  che  k0   Flint,  1  hey 

obtained  by  computing  the  wel     ii    a|  ince  of  the  </   rlo 
components,        '.•■/,    the    aircraft  was  wel  jhed  In  various  coi  flgur  - 
;  wit'  the  fusi  I  age  reference  Ii  i  horizontal  ,     rhree  p  I  ji  f 
iles  placed  under  the  I  indi  ig  gear.  This  method       r  is  weight 
and  horizontal  e.g.  position  o\   ly.   Thirdly,  the  alrcr  ft  was  h 
froi  'J  beam  by  dynarra  m   r  an  i  c   I-   itt  iched  to  the  rotor 

shafto   Foi  various  confj  uratlons  fhe  angle  ot  Incll   tlon  of  the 
fuselage  ref  rence  line  and  weight  were  determined.  Knowing  ttie  In- 
clination of  1  he  rotor  head,  a  check  col  Id  be  ob1  lined  oi      vertical 
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fa  I  :«  .position.  The  e.g.  I  oca  1 1  n  fh   obtained  i 
;   i  a   rate  to  with  in  j£  0.25  inch  Ii  l"he  horizon  to  I  Hi   f  ioi 
±  0.50  1     i   f he  vert  I ca I  direct  ion.   Th   .g.  pos It  ion  (foi 

eighl  tested)  re  I  itlve  to  ("he  rotor  head  Is  Ii  ,ted  ii       I   I  . 
It  m      well  to  nolo  that  since  the  aircraft  [s  relatively  light, 
Co  j.  pos  M  Ion  can  be  affected  •  ,    hi  w  the  pi  lo1  sit  ;  In  the    it.  .  urli 
Plight  festing  every  efforl  was  made  to  malntai   the  c.  .  po  It  ion 
determined  above  and  ,ross  weight  is  felt  to  bo  accurate  to  within 
+  3  poui  j  . 

Rotor  Blade  Section  Characteri  ii 

To  accui  i   ly  determine  the  section  aerodynamic  propi  i  i    jf  1  he 
rot<  i-  blade  a  12  ircii  section  of  the  blade  was  te  ted  In  ("he  J  im< 
Forrestal  Research  Center  tw<  -  I]  i  nsional  wind  tunnel.  Th  re   I  s  of 
the  wind  funnel  test  are  plotted  or;  Figure  19.   Tfn-  section  Irng  oeffl- 
cient  was  corrected  for  Reynold;  i    i  i  el      to  the  average  operating 
condi  i  ns  a1  fhe  0o7M   a  d  expressed  as  a  quadratic  i:     i'l      I   of 
attack.     .  jot:  lift  curve  slope,  lr   coefficient,  and   I  ide  pltc 
ie  thus  obtaj  led  are  rabul  r  I  in  Appendix  I. 

Proce J jre 

The  f  I  i  gh1  testing  procedure  followed  was   imply  To  e  i  il  llsh 

steady-state  level  flight  conditions  at  five  knot  intervals  from  30  to 

0:0  knots  ii     fed  forward  speed,  i  d  record  indicated  spei  ,  r  for  r.po-  ., 

control  arigle  and  fuselage  angle  of  attack.  Since  J  I  festlng  was  lone 

at  altitudes  of  from  200  to  300  feei  M>L,  no  tei  |  asured  i 

flight  and  density  corrections  were  applied  using  corrected  burface 
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rhe     olti     i  l  is   vo  1 1  >■  i    v  curve  thus  obt.i ii 

is   plotted  on   Fi  jure   20,    fh     rotoi     refers  i     ■  lie   of     it1  ick 

velocity  on    FI jure   2\ ,    and   angle   of   attack   of    the   no    feathering   pine 
versus   velocity  on   Fi  ;ure  22.     The   expei  i      ntal    hhi       I    coefficient 
versus    iip    speed   ratio  obtained    is    plotted   on    FJgure   8    for   compri    oi 
witli   theoretical.      Each   point   plotted   represents   an   averjge  of    at    |i 
two  readings. 
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IV.   REJULT.l  ,ND  Dl  CU  SI  ON 

The  re  hor  foi  e  and  moment  derivati  .     /e  >eei   ij  ul  -tod 
led  condi  t  ions  of  constant  rotor   ;  r  I  able  rotor 

speed  in  a  -       .  i.h  torque  equjll  ri  im,  and  are  presented  in 
Figures  6,  1,    10,  II,  13  ai  j  |4.   The  conditions  of  constant  r^icr 
speed  and  vai  i  il  le  rotor  ;peed  in  accordance  with  torque  equilibrium 
are  fell  to  be  limii  |ng  c  >i  li  i  is  f<  r  jut     I  1  ve  rotor  operatlono 
The  assumpl  i  i   I  con  t  n  t  rotoi  :  poo  !  ould  be  [i  herpreted  to  m 

the  Inert  1  ■  of  I"  he  rv  tor  is  verv  I  r  j'  (although  it  i   usually 
justified  by      |ng  that  rotor  poo  J  cha  ges  orv   negllgi  ly  srrull); 
\  i  i  .  le  re  I  ;   pe  j  was  obtained  through  the  assumption  of  negllgli  le 
i  erti  io  Obvjou  ly  l"he  real  rotor  has  sumo  Inertia  a  .d  its  operation 
will  lie  somewhere  between  the  two  c^ses.   For  slow  perturl  itloi 
rotor  should  hei  d  to  act  more  like  I       IJglble  inertia  cose;  for 
rapid   changes  inertia  will  tend  to  maintain  a  constant  rotor  speed. 
The  question  -  i  irally  arises  then  js  to  how  rapidly  the  rotor  accelerates^ 

Isky  ii-  i  ites  ii  w  rk  with  i  helicopter  rotor  of  ;imj  I  ir  -olldity  In 
autorol  i1  i    (  v  f.   )  that  the  re  "or  i      <    Its  peak  acceji  r  it  Ion  In 
about  2s  seconds.  Using  ihis  figure  on  \\  >  .     er  of  r.   I    .it 

may  he  concluded  that  tor  long  period  dynamics  the  assumption  of 
negligible  Inertia  for  rhe  autorotatlng  rotor  is  v j } Id. 

In  the  longitudinal  p  I  .j  e   I  hie  stability  Is  usually  leflned 
as  the  restoring  moment  In  relation  to  a  change  In  forward  linear 
velocity  at  constant  angle  of  attack  and  In  relation  to  a  change  in 


irigh   i         I    nstant  ve  \o>   Ity.  A  comparison  of       loclty 
I  ibllity,  NL,  (Figure  I  )  shows  that  the  in  rtl  ile  ;  roti  r  i   ipproxl- 
>ly  neutrally  stable  while  l"he  constanl  ;peed  rotor-  ]-■   quit'   I   I  „ 
A  comparj:  ngle  of  attack  stability,  M.  (Figure  \  \) ,    shows 

that  the  autororatlng  rotor  becomes  strongly   I  id le  iue  t< 
changes  whi  le  the  constant  speed  rotor  Is  unsl  ib|e.   The  st  itl     ibilltv 
derivatives  of  1  he  autorotatlng  rotor  of  low  Inertia  then  behave  r  >ore 
I  Ike  a  fixed  wing  than  a  helicopter  rotor « 

To  qualitatively  explain  how  l~he  ro\    r  speed  vari  .lion  affe< 
the  sialic  stability  consider  a  rotor  ifi  trim  forward  fii  hi.  An  In- 
crease in  forward  velocity  with  constant  ai  jl<  of  attack  will  c  use  > 
sroa  I  I  Increase  in  the  rotor  resultant  f  jr:e.   Further,  the  Increase  In 
velocity  will  increase  the  longitudinal  fl  ippli  ;  angle  causing  i  nose  ip 
or  positive  moment  from  trim  If  the  rotor  speed  is  constant  a   The  rotatl*   I 
speed  of  the  autorotatlng  rotor  Increases  with  velocity;  however,  which 
causes  the  longitudinal  flupphq  angle  to  decrease,  resultii  j  in  a  r  ose 

or  negative  moment.  When  the  rotor  speed  variation  with  forward 
speed  occurs  so  that  the  equilibrium  forque  Is  maintained  the  net  result 
is  10  variation  in  longitudinal  f|  ipping  angle,  and  no  m    I  is  pre 
except  due  to  1  he  small  increase  |n  rotor  resultant  force.   The  >uto- 
rotatlng  rotor  is  therefore  approximately  neutrally  stal  le  with  respe  I 
to  speed.  This  fact  could  have  been  deduced  by  examining  the  aero- 
namlc  torque  relation.   Since  the  torque  can  be  expressed  as  a 
tlon  of  u,  and  a,  It  could  further  be  expressed  as  a  function  of 
p,  and  a  by  appropriate  substitution  from  momentum  ind   I  ide  e|em< 
theory  for  \  and  Ct.   Since  by  assumption  the  equJII  rium  Torque  I 
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maintained,  it  then  follows  that  the  tip  speed  ratio  i       ■  I"  1 

pei   i   it  ions  a1  constant  angle  of  attick.  And  since  the  longitudinal 

flapp]     gle  can  also  be  completely  expnessed  as  a  fun<  ion  of  p,  and 

y  i he  same  s  istjtutlons,  it  too  must  remain  constant  for  pertur:  i- 
tions  at  constant  jngle  of  attack*  The  method  for  quant it-.it ive I y 
pre.ju.i  ig  the  velocity  stability  of  the  autogyro  developed  in  the 

alysis  section  did  not  yield  Identically  neutral  velocity  stability. 
Since  all  1  he  mome'  ts  acting  are  aerodynamic  (power  effects  have  been 
no    ted  in  velocity  perturbations)  then  all  the  changes  In  moment  at 
constant  angle  of  a_i  tack  are    functions  of  velocity  and  the  net  result 
should  have  been  identically  neutral  velocity  stability*   It  is  felt  that 
the  assumption  of  negligible  fusel  ige  moment  and  an  overest Imat Ion  of  the 
tail  contribution  accounted  for  this  apparent  error. 

To  see  the  effect  of  variation  of  rotor  spec  !  i     |le  of  attack 
stability  consider  again  the  rotor  In  trim  forward  flight.,   Increasing 

jle  of  attack  and  maintaining  constant  forward  velocity  will  cause  an 
increase  in  the  rotor  resultant  force.  The  longitudinal  flapping  angle 
will  also  increase  (for  constant  r.p.fru)  tending  to  produce  a  further 
nose  ip  or  unstable  moment.  However,  the  increase  In  angle  of  attack 
will  cause  the  autorotiting  rotor  to  Increase  rotor  speed,  resulting  In 
a  decrease  in  the  longitudinal  flapping  angle  and  a  further  increa.e  i 
resultanl  force.  The  net  result  for  the  autorotatlng  rotor  is  to  produce 
a  nose  down  or  stable  moment  due  to  the  strong  tendency  for  rotor  speed 
to  increase  with  angle  of  attack. 

If  the  coefficient  of  the  constant  term  of  the  characteristic  or 
,  ii    "  equation  is  used  as  a  measure  of  static  stal  I  I  Ity  (i  neces  .  r, 
it  not  suffi<  ion.  condition  for  dynamic  stability): 
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E  =  Wcos-YCZyM^-Z^jMy)    C  for  si  it  I  c  st  i  II  Ity 

It  Is  observed  :  hat  the  aircraft  Is  statically  the  assumed 

conditions  of  both  constai  I  ind  var  ia  ►  I e  rotor  ;peed.  As  hown  on 
Figure  23  the  static  stability  of  the  alrcr  iff  is  I  iproved    v.ri   I 
rotor  speed  due  to  the  Improvemenl  |i  i  hoi   in  ||e  of   ttacl   I   Mi.  y. 

I  t  can  be  '.how:i  th  -it  hh<   hat  Ic  s  l"ab  I  I  I  ty  of  hhe  a  I  re  raft  i 
Independent  of  horizontal  e.g.  position  provided  thai  [bo    .  .>  travel 
Is  limited  to  the  range  where  adequate  control  I   ivalla  |e. 

.t  I  tut  3  ig  equations  (15a)  and  (15c)  Into  the  static  it   I  I  Ity, 

or  E,  tern  gives 

cWfp     .Mfp 
E  =  Wcosy  ''    z(Z  Hy-EyH  )  +  Zy  ~fa  Z0  ~       ' 

and  is  independenr  of  x.   This  could  be  ph   I   lly  reasoned  by  not  1 1 
the  trim  control  angle  6  is  a  strong  function  of  the  hoi  izor;j|  c.ge 
positic;-  through  the  pitching  moment  equlil  rlum  equation,  equatloi  (40) , 
and  hence  vary  in   h  rizonta]  e.g.  position  will  vary  the  trim  control 
angle,  bu1  The  static  stal  illty  will  not  change,  Varyln  |  The  horizontal 
c8g.  position  will,  of  course,  affect"  th   I  italllty  contributions  of 
fuselage  and  power  by  varyli  g  their  moment  arms  but  this  i  :.  a  secondary 
effecto   Static  sta  ilir/  can  :e  improved     ovlng  the  c.  .  vertically 
but  is  physically  impractical  foe  an\   I  nlficant  Improvement* 

Trim  curves  of  control  at  jk,  attitude  angle  ,    a      I   of 
attack  of  the  no-feathering  plane  respectively  versus  forward  velocity 
taken  from  flight  testing  are   presented  on  Figures  20,  21  and  22.  A 
comp.jr  i  son  of  their  measured  and  predicted  slopes  at  various  forward 
velocities  is  presented  In  Table  I  he  low. 


Tab|«     ! 
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V    f I    per    sec 


A' 
AV 

Act 
AV 


AV 


60 


Predicted   '  -0.000C77 
Fl Ight  test  :  -0.0007oo 


70 

-0.000C.44 
-0.000684.; 


Predicted    -0.00731  »  -0.00;0| 

Fl  1  ;!  t  tost  ;  -0.00833  |  -0.00 567 

Predicted  |  -0. 00849  -0.00 

Flight  test  j  -0.0085  j  -0.00o4:> 


80         90 

-0.000640  -0,000648  ■  -0.00C 
-0.000617  ;  -0.000434  '  -0.000435 

-0.00344   -0.00248   -0.00|79 
-0.0035   i  -0.00234   -0.00142 

-0.00403  !  -0.00  I    -0.00244 
-0.00425  i  -0.002' -7  |  -0.001  i 


Considering  the  rather  unsophisticated  type  of  Instrumentation 
used  i  ri  flight  testing  and  approx  Imat  ions  necessary  In  the  analysl  ,,    I 
correlation  between  predicted  and  measured  data  appears  good.   It  i 
much  better  at  lower  forward  velocities  than  high,  duo  pro  ably  to  i 
weaker  Influence  of  the  unknown  and  neglected  fuselage  moments. 

The  perturbation  equations  of  motion  give  us  a  relation  between 
flight  variables,  such  as  V,  a,  and  6,  and  the  stability  derivatives  whl  h 
describe  how  the  forces  and  moments  act  on  the  aircraft  and  which  wo  wl  h 
to  determine.   Since  each  independent  relation  necessarily  contains  m,:y 
of  the  fli  jht  variables  It  Is  seldom  possible  to  obtain  an  explicit  re  I  i- 
tlon  for  any  one  unknown  stability  derivative  In  terms  of  the  flljht 
.  irlables.   It  Is  possible  In  wind  tunnel  testing  to  hold  all  flight 
variables  but  one  constant  and  measure  the  resultant  change  In  for. 
or  moments  that  act.   Flight  testing  methods  cannot  hold  all  fli 
variables  but  one  constant  and  measure  the  resulting  changes  because 
the  aircraft  cannot  physically  satisfy  the  equations  of  motion  under 
these  conditions.   As  a  result  steady  state  flight  test  methods  usually 
attempt  to  hold  as  many  flight  variables  as  possible  constant  and  verify 


i  I  i  I  y  j'  i  i  •/  ,  i  .  i    hi     i   los  of  thi     r  I  .  i   .   Tl  i 
in  I  [s  analysis  by  writing  the  pertur  it  Ion  iquatlon1  I 
Imple  a  form  as  possible  and  still  retaining  i;    i  nlflcant  derlv  tl 
Tl.     ult  Is  1  he  two  co.-trol-f  i  xed       -  I    pertui   1  Ion  equ  .tic 
(I  a)  md  (16b) 

V  +  :aAd|  +  Z^y    m   0 
My  V  +  M.A-|  +  M6:3  =  0 

All  the  ma  joe  stability  derivatives  (with  the  exception  of  pitch 
damping,  a  rate  derivative)  are  contained  In  these  two  relatlonSe   By 
ellml   -I  ,  any  one  of  the  flight  variables  between  the  two  equations 
a  relation  Is  ol  Mined  for  the  stability  derivatives  In  terms  of  ratios 
of  the  two  remaining  flight  varjables«  The  ratios  predicted  In  Table  I 
were  obtained  from  the  above  relations  as 

Z^MV-ZVM^ 


IV   Za(Ma-Mfi) 


or 


Ach    .'vM'"-.:-Mv 
AV    Za(Ma-M6) 

A  further  check  on  the  2  force  derivatives  Is  obtained  from  the  Z  force 
equation,  noting  that  *a  =  Aa,  +  A5 

Act    Zv 
AV  "  Zc 


M  -'•'  had  known  one  of  the  Z  force  derivatives  fron      [i  lepci 
source,  su  h  as  a  wind  tunnel,  then  the  ether  could      ptl  Itly 
verified.  Although  steady-state  f  1 1  gh1  I       hnlques  have  thl  i  .;  I  .- 

,  ("heir  ad  van  1     of  dl  associatI<  n  with  phase  lag  prob  lei 
and  use  of  simplified  In  trumei  I  itlon  still  make  them  attract  I  vi  In 
many  cases.  Often,  by  selecting  proper  if  I  hechnl  ;!i  .,  |1  Is  possible 
to  cross  plot  data  to  evaluate  derivatives  not  otherwise  obtalm  I  I  , 
as  Is  often  done  in  fixed  wing  aire  rift  for-  the  pitching  mi  ,    -  i  jle 
f   1  ick  stabjllty  derivative. 

The  flight  test  plot  of  control  angle-velocity  curve  (Figure  20) 
shows  a  gradual  decrease  In  the  negative  slope  with  forward  velocity. 
The  predicted  slope  of  this  curve  is  essentially  consl  int.   It  Is  feli 
that  the  influence  of  the  fuselage  momer.ls  causes  this  decrease  i 
slope  aT  higher  velocities,  and  neglecting  these  moments  |i  M      lysl 

resulted  In  the  Inability  to  predict  the  change.  The-  concept  of 

t 
"stick-fixed"  stability  as  used  In  fixed  wing  airplanes  implies  a 

definite  relation  between  stick  position  and  forward  speed;  that  is, 

the  stick  position  corresponding  to  a  higher  forward  speed  m  ;i 

farther  forward  to  be  regarded  as  stable.  A I  thou jh  the  corollary  Is 

not  necessarily  warranted  for  rol  iry-wlnged  aircraft,  the  autogyro 

tested  Is  stable  by  this  criteria.   The  fact  thai  the  slope  decreases 

rather  than  Increases  with  velocity  could  indicate  a  possl  i|e  Instability 

at  some  high  forward  speed.  Hufton  et«>a|,,  in  an  Investigation  of  a 

British  direct-control  autogyro  (Ref.  8)  reported  suc-h  a  high  speed 

Instability,  but  determined  the  cause  as  blade  twist. 
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The  altitude  angle-velocity  and  angle  of  attack  of  Ihe  r,o— 
feathering  plane-velocity  curves  also  show    I  i  correl  it  I  n  at  the 
lower  forward  speeds,  but  an  estimation  of  the  accuracy  of  verifying 
■Ihe  stability  derivatives  is  difficult.   |n  all  cases  below  90  f.p.s. 
the  predicted  slopes  were  within  ]2%   of  the  f  I  J-ght  test  curves,  which 
themselves  undoubtedly  contain  some  measurement  error*.     The  ritlo  of 
Zy  to  Zq,  verified  by  the  slope  of  the  angle  of  attack-velocity  curve 
would  easily  point  up  any  significant  err^r  in  either  one  of  the  deriv  - 
tlves,  ;.  ut  if  the  error  existed  in  both  derivatives  by  a  constant  fjctor, 
1hen  the  error  could  go  undetected..   Significant  discrepancies  In  the  pre- 
dictions of  the  control  angle-velocity  or  altitude  angle-velocity  ratios 
would  be  more  difficult  to  analyze  Assuming  the  Zy/Z  ,   rjtio  was  verified, 
each  ratio  still  contains  the  three  moment  derivatives..  Due  to  the  re  I  itlve 
magnitudes  of  these  three  derivatives  (My  Is  very  small  compared  to  Ma  and 
Mf.)  a  large  error  or  even  a  change  in  sign  could  exist  in  My  and  would  be 
difficult  to  verify.   In  fact,  assuming  neutral  velocity  stability  for  the 
aircraft  as  previously  discussed  improves  the  correlation..   Since  the  M$ 
derivative  is  essentially  the  rotor  M(y  plus  a  constant  term,  the  control 
angle  and  altitude  angle  curves  essentially  verify  the  pitching  moment- 
angle  of  attack  stability  derivative,  which,  as  pointed  out,  becomes 

increasingly  difficult  to  predict  at  hljher  velocities. 

V 
In  summary  it  is  felt  that  the  trim  curves  tend  to  verify  the 

vertical  force  derivatives  and  angle  of  attack  stability,  and  can  point 

out  possible  unstable  areas  not  predicted  by  analysj  o   Neutral  velocity 

stability  for  the  aircraft  is  felt  to  be  a  reasonaole  approximation. 
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V.   COf.CLU.ION, 

As  a  resiiM  of  this  analysis  the  follow!.,  ;  conclusions  art' 
f  e  I  f  to  De  va  I  I  d : 

1.   The  method  developed  in  this  analysis  for'  predict) 

static  stability  derivatives  is  valid,  but  is  limited  y 

fuselage  effects. 
,   2.   The  autorotatinj  rotor  i  s  approximate  Jy  neutrally  stable 

with  respect  to  velocity  chanjes. 

3.  The  autorotating  rotor  is  statically  stable  with  respect  to 

altitude  change. 

4»  The  static  stability  of  the  direct  control  typo  ro ror  is 

Independent  of  horizontal  e.g.  position  as  long  as  it  is  wlthl 

contro I  I Imlts. 

5o  The  direc"!  control  type  autogyro  considered  is  statically 

stable,  due  in  large  part  to  the  a  ility  of  the  rotor  to  change 

speed  during  pertur  at  ions. 

6.  Trim  curves  tend  to  verify  ("he  method  developed  for  conput  i 

static  stability  derivatives,   j    /  themselves  are  not  conclusive. 
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PEN   I  >     I 
», 

P ,y  .  1     ;  1   Ccn  •  yi  o 


Gross    w.  i  )h  r 


N    =  41  C    I    . 


C-  ntei    of     iravll  y    In    rel  it  Ion   to 

i     tor    reference    axi  s 

f 

Equivalent  flat  plate  dr    rea, 
fusi      (est  ]  nated) 

Engine  -  4  cyl .  M  Cu I lou  h  4318-A 

Thrust  axl   p  rallel  to  fuse  la 

reference  axis  (a  re.,) 

P      f 

Propel  ler  -  45-24  Ba  ik   -Maxwel  I    two-:  laded   pu  ,her 

Hor  I  zonl  1 1  hai  I  -  Ai 

i"  1  o 
Lift  curve  s  lop- 
Sweep  of  quarter  chord 
Till  arm 
T  :  I  I  i  nc I dence 

Roi or  (teeter  hinge) 

[I  ;  i  er  of  •>  I  ados 
Radius  of  b I 
I  I  ide  chord  (con  tant) 
Rot  »r  so  I  idlty 

Bui  |t-Ii in  ;  ang  le 

I !  I  ide  p  1  tch  angle 

i  i  ide  twi  st 

;  i I  jde  airfol I  section 

Li  ft  curve  s lope 

Drag  coefflclei   is  a  power  ;erl 


Mean  dr  ig  coe f f Ic lent 
M    moment  of  Inert  I  i  about  the  fl  pplnj 

h  i  pge  (pc (   I  . Je) 
Mass  moment  about  f  I  jpp I rg  hinge 

(per  blade) 
Weight  moment  about  flapping  hi 


;acR' 


(per    Made) 
I  '  s    b|  ide   constant   y .    - 

■  I 
Incidence  of    plane   perpend  I c j I  ir    to   rotor 

reference   axis   to    fuselage    reference    line 


x   =  -0. 10      ft. 
z    -  -    .1    5    ft. 


f   =  8    ft 

72  IP  (  I) 

xp    =  -    .         f  i 

Zp    =  +0.004|o    ft. 


S+    «  2. 
A      =    I .54 

=  2.0 
..   =  4  P 

--I.82    f 


t 


=  0 


b   = 

R    m     |0    ft 

:    ■  0  o :    2    f  1 
c  «  G.03^3 
;c=  0.07   i   idlans 
eo=  0.02-     r    Hans 
9.-  0 
B<  t,  ,enG~3 
a   =   ■  o  j  per   radian 

=     c  +    |-r  +  2>2 
cc   m  CO  I 
■ .    -  -0.0509 

\   =  0.3170 

"   =  0.01 

I  j  =  20.77  slug  ft' 
M  =  2.87  slug  ft2 

M,.;  =  92.   lb  ft 

■   =  i.  A 
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APPENDIX    EL 
Aerc  lyn   mlc    Tor  ;  I     i  ion' 


Bailey,    in   Ref.    4,    jives   the-    following   rejotions    for   the 
accelei    itingand   decolor   ting   torque   actln      o  lifting   rotor    h<w| 

no  fwlst:      (neglecting   terms  of   the  order  u/    and   higher) 


I.      Acce  lerat  in>j  Torque: 

2CQ       ,  2  : 

—  -  ^t^,)^+a(t^2)XG0+a(t4/4)% 


where  B2         U       Y)2B8\ 

B3     8B   y  2B9\  ? 

ty|   O  ■  +     +    1   .    u/ 


4>2     '   3      5     864 


8B2   v2B10   „ 

t    =       +  J .  u- 

4,4  1080  I 


2c  Decelerating  Torque: 


-^-  6o<t5,i)+6ir  <+5,2>*+<+5,3>Qo  '^2r  Ct5^>x2+(t5^)  e0+<t3>8)e02  ] 

where  ,  , 

2 

f5,|    "  4      +  4  H 

t  -I 

5*2        3 

< 

1         '      2 
f5,3  =  4"  +  4   M- 


'I  II  I  2      /B4  BD        B6 

t5.5=-+     --+^+-T  +  Y.         , +  n 

'2  4        L  2B4  '  162        81         144/ 


n 


2  4  4  B-  Bl         B7  I 


,   =T  + 


,  a  tip  loss  f jctor  B  «  I-  —  =  0.972,  substjtul ;      ?  known  rotor 
erlsl  ics  as  listed  in  Appendix  lit'   hhe  ibove  equations,  jnd 
addi  lerating  to  decelerating  torque  yield  thi  following 

expi   .si  in  for  the  total  aerodynamic  torque: 


Cg   =-  r    m.2(AjX2-M2A.4A3)    -    (B|X.2-te2^-fe3)    n 


v.h<.  re 


A,    =      -  .47  1 
A2    =      C.  5jCA 

=  -0.001 144 
3,  =  -2.429250 
B2  =  -0.035330 
B5   =     0.004121 

At   an   equlll    rlum   rotor-    speed,    the  torque    is   zero, 


CQ   =  0 


2  B|\24B2'\-+e3 


»    "    TT 


A  in    ^A-fnj 


F^r  J   i  /e      iip  speed  ratio  the  above  relation  results  In  a  quadratic 
ii  Inflow  ratio.   The  I  r  }er  arithmetic  value  obtained  from  the  solution 
of  the  quadratic  corresponds  to  rotor  operatloi   1  -  positive  am  I 
of  attack. 
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